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We review some recent developments in the theory of rotating atomic gases. These 
studies have thrown light on the process of nucleation of vortices in regimes where 
mean-field methods are inadequate. In our review we shall describe and compare 
quantum vortex nucleation of a dilute ultracold bosonic gas trapped in three dif- 
ferent configurations: a one-dimensional ring lattice, a one-dimensional ring super- 
lattice and a two-dimensional asymmetric harmonic trap. In all of them there is a 
critical rotation frequency, at which the particles in the ground state exhibit strong 
quantum correlations. However, the entanglement properties vary significantly from 
case to case. We explain these differences by characterizing the intermediate states 
that participate in the vortex nucleation process. Finally, we show that noise corre- 
lations are sensitive to these differences. These new studies have, therefore, shown 
how novel quantum states may be produced and probed in future experiments with 
rotating neutral atom systems. 
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1. Introduction 



Vortex nucleation is a topic at the heart of the nature of superfluids and their 
intrinsic quantum character. Superfluid flow is a direct manifestation of quantum 
mechanics at the macroscopic level and is only stable below a critical velocity. 
Above the critical velocity the generation of phonons, rotons (in 4 He) or vortices 



can lead to a breakdown of superfluidity ( |Ihas et al. 1992). The microscopic nature 
of these processes cannot, however, be studied in superfluid liquid helium. There 
has been a resurgence in the study of vortices and their production in the new 
class of superfluids produced using ultracold atoms. Indeed ultracold atoms offer 
a unique opportunity to investigate topological excitations as recent experiments 



have demonstrated both with bosonic and fermionic atoms (Madison et al. 2000 



Raman et al. 2001 Haljan et al. 2001 Zwierlein et al. 2005). The production of 



vortex arrays has, in fact, been crucial in demonstrating presence of a superfluid 
order parameter in the ultracold atom systems. 

Aspects of vortex nucleation in ultracold atoms have been analyzed in various 



theoretical studies based on a mean-field treatment (Feder et al. 1999 Sinha & 



Castin 2001 Kasamatsu et al. 2003 Lobo et al. 2004). More recently, attempts 



have been made to approach situations where the neutral atom vortices require 
methods beyond mean-field techniques. This should be the case, for example, in 
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rapidly rotating clouds where the neutral atom vortices are analogous to those 
studied in the quantum Hall effect. An overview of the experimental field has been 



given by Stock et al. (2005), while Cooper (2009) reviews the different regimes 



theoretically 

In recent studies, the strongly-correlated nature of the ground state around 
the region for quantum nucleation has been elucidated. The importance of such 
studies is at least twofold: firstly production of novel quantum states and secondly 
investigating the role of strong correlation in the nucleation process. Two different 
physical scenarios were considered: In one of them the superfluid ultracold atoms 
were trapped in an optical ring lattice that can be rotated to introduce angular mo- 



mentum into the system 1 
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; in the other the atoms were confined in 



a rotating asymmetric trap under conditions equivalent to having charged parti- 
cles in a magnetic field ( Dagnino et al. 20096|q I . In the case of the rotating ring 
lattice the authors were most concerned with the production of Schrodinger cat 
states, i.e. macroscopic superpositions of states with and without a vortex. They 
were inspired by the similarities of this system with a superconducting quantum 
interference device (SQUID) which exhibits macroscopic tunneling between states 
of opposite current flow (Rouse et al. 1995). As discussed in Leggett (2002), these 



macroscopic superposition states are important for testing the limits of validity 
of quantum mechanics and can be used to achieve quantum-limited measurements 



in precision spectroscopy (Leibfried et al. 2004 Cappellaro et al. 2005 Lee & 



|Khitrin| |2005 ) which is important for ultra-precise gyroscopes. In Hallwood et al. 
(|2006[ |2007j )7 |Rey et al\ ( |2007[ ), |Nunnenkamp et al.\ ( |2008| , and |Nunnenkamp k 
|Rey| ( |2008| ) the goal was to find the optimal conditions for cat state production. In 
the case of the vortex nucleation in the lowest Landau level (LLL) Dagnino et al. 
(2009& aj demonstrated that the mean- field picture breaks down close to the nucle- 
ation point. The nature of the correlated states, while highly entangled, was clearly 
distinct from those produced in the rotating ring lattice. 

In the present work we show how this difference comes about by studying the 
nature of the intermediate states, i.e. the routes that connect states with and with- 
out vortices. By constructing effective Hamiltonians, we clarify the competing role 
of interactions and trap or lattice asymmetry in each of these cases and show that 
the entangled nature of the ground state and the strong correlations involved in the 
nucleation process can be seen in quantum noise correlations. We believe that ro- 
tating atomic gases offer significant new opportunities to study strongly-correlated 
atomic systems. They may also have applications to areas where such states can be 
used in quantum information science and precision measurement. 



2. Vortex nucleation in the rotating ring lattice 

In this section, we investigate the effects of rotation on ultracold bosons confined 
to one-dimensional ring lattices and superlattices. This is an attractive system for 
study as it can be produced in the laboratory and isolates important issues of 
the underlying physics of the nucleation process. We find that at commensurate 
filling there exists a critical rotation frequency, at which the ground state of the 
weakly-interacting gas is fragmented into a macroscopic superposition of different 



quasi- momentum states (Hallwood et al. 2006 2007 Rey et al. 2007 Nunnenkamp 
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et al.\ 2008 Nunnenkamp & Reyl 120081). We note that iWatanabe & Pethickl (]2007h 



and Danshita & Polkovnikov ( 2009 ) have studied related aspects in a similar system. 



(a) Hamiltonian 

We consider a system of N ultracold bosons with mass M confined in a ID 
ring lattice of L sites with lattice constant d. Optical ring lattices can be created 



with Laguerre-Gaussian (LG) laser beams, as proposed by Amico et al. (20051 and 
realized by Franke- Arnold et al. (2007). LG beams can be derived from ordinary 



Gaussian beams, e.g. by means of computer-generated phase holograms (Chavez 



Cerda et al. 2002 ) . A conceptionally different approach to arbitrary 2D potentials 



are spatial light modulators which have recently been used for dynamical manipu- 
lation of Bose-Einstein condensates (BECs) ( |Boyer et al. 2006 ). The ring is rotated 
in its plane (about the z-axis) with angular velocity fi. By subtracting the rotation 
energy H to t — J dx &VLL Z § we transform to the rotating frame where the potential 



y(x) is time-independent and the many-body Hamiltonian is given by (Bhat et al. 
2006||Rey efa^|2007" ) 



H = / dx$ t 



2M 



V(x) 



4nh 2 a s 
2M 



(2.1) 



In this expression, a s is the s-wave scattering length, V(x.) the lattice potential, L z 
the angular momentum and x the 3D spatial coordinate vector. <3E>(x)^ and $(x) 
arc bosonic creation and annihilation field operators. 

We assume that the lattice potential V(x) confines the motion along the z-axis 
as well as the radial direction in the x — y plane so strongly that only the motion of 
the atoms along the ring has to be taken into account. In addition, we assume that 
the lattice is deep enough to restrict tunneling to nearest-neighbor sites and that 
the band gap is larger than the rotational energy. These assumptions imply that 
the bosonic field operator $ ca n be expanded in W annier orbitals confined to the 
first band $(x) = £\ ajWj{x) (Jaksch et al. 1998). Here, Wj(x) are the Wannier 
orbitals in the rotating frame and aj the bosonic annihilation operator of a particle 
at site j. We recall that the Hamiltonian of a neutral particle in a frame rotating at 
frequency fl around the z-axis, H = p 2 /2M — QL Zl is equivalent to the Hamiltonian 
of a charged particle in a magnetic field along the z-axis, H = (p — A) 2 /2M with 
the effective vector potential A(x) = MQ,(z x x). This implies that we can first 
calculate the Wannier orbitals of the stationary lattice Wj (x) and then account for 
the presence of the effective vector potential A(x) via the gauge transformation 

Wj(x) = exp /* A(x') • dx'] W,(x). 

In terms of these quantities and up to on-site diagonal terms, the many-body 



Hamiltonian can be written as ( Bhat et al. . 2006 Rey et al. 2007 1 



3 = 1 



In this expression, fij = a^aj is the number operator at site j, 9 is the effective 
phase twist (Peierls 1933) induced by the gauge field, 8 = f* i+1 A(x') ■ dx! — 



(2.2) 



3=1 
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, Jj is the hopping energy between nearest-neighbor sites j and j + 1: 



'h = I dxW* 



2M 



V 2 + V(x) Wj+i, and U the on-site interaction energy: U 



(b) Superposition states in rotating ring lattices 

We start our discussion with the case of a uniform ring lattice, i.e. Jj = J for 
all j. To understand the effec t of r otation on the atoms in the ring lattice, we write 
the many-body Hamiltonian (2.2) in terms of the quasi- momentum operators b q = 
,-2mqj/L ^ wnere 2nq/dL is the quasi-momentum and q = 0, . . . , L — 1 



j_ 

sft 



an integer. In this basis the Hamiltonian (2.2 1 has the form (Hallwood et al. 2006 



Rey et al, 20071 



L-1 

H = H sp + H int = E E q b\b q 

9=0 



U 
2L 



L-1 

"'/"•- 

q,s,l=0 



blbthb^ 



i-i] i 



3dL 



(2.3) 



where E q = — 2J cos(2Trq/ L — 8) are the single-particle energies, and the modu- 
lus is taken because in collision processes the quasi-momentum is conserved up 
to an integer multiple of the reciprocal lattice vector 27r/d, i.e. modulo Umklapp 
processes. 

Following |Hallwood et al.\ ( |2006[ ) and |Rey et al.\ ( |2007[ ), we show in Fig. [l] the 
single-particle spectrum as a function of the phase twist 6. In the absence of rotation, 
i.e. 6 = 0, the state with zero quasi-momentum \q = 0) is the single-particle ground 
state. In the rotating system the ground state depends on the phase twist 9. Writing 
6 = ^-m + ^ with m an integer and < A9 < 2n, the ground state is the quasi- 
momentum state \q = m) for < A9 < it and \q = m + 1) for tt < A9 < 2ir. 

In the absence of interactions U = 0, the ground state of a bosonic many- 
body system is the state with all N bosons occupying the lowest-energy single- 
particle state. At rotation frequencies corresponding to phase twists with A8 = tt, 
the ground state of the single-particle Hamiltonian H sp is two-fold degenerate, 
i.e. E m = E m+ i, so that there is a N + 1-dimensional degenerate subspace at the 
iV-particle level (see Fig. |2]). A convenient basis for this subspace are the Fock 
states |n, N — n) with < n < N, where n particles are in the quasi-momentum 
state \q = m) with energy E m and N — n particles in the quasi-momentum state 
\q = m + 1) state with energy E m+ i, respectively. 

For weak interactions, i.e. U <C J , we can use first-order perturbation theory 
to account for the effect of interactions. It predicts that the energies of the states 
In, N — n) are 



U 



= nE m + (N- n)E m+l + — [A^ - 1) + 2n(N - n)] 



2L 



(2.4) 



We see from this expression that the degeneracy is lifted and the states of lowest 
energy are \N, 0) and |0, N). These two states are still degenerate and higher-order 
coupling is needed to break the degeneracy (see Fig. [3]) . 

At this point it is important whether the number of atoms N is commensurate or 
incommensurate with the number of lattice sites L. While there are many different 
paths that couple the states \N, 0) and |0, N) in the commensurate case, in the 



Article submitted to Royal Society 



Review. Quantum Vortex Nucleation. 



5 



incommensurate case there is no coupling between these two states and they thus 
remain degenerate at all orders of perturbation theory. To understand this fact, let 



us consider the total quasi-momentum operator K 



■2ts_ 
L 



Since 



mod L 



the many-body Hamiltonian ( 2.3 ) commutes with the quasi-momentum operator K, 
i.e. H, K = 0, the Hamiltonian has block diagonal form if the quasi-momentum 

Fock states are ordered according to the eigenvalues of K. In the commensurate 
case, we have N — nL with the number density n being an integer, so that the 
states \N, 0) and \0,N) have total quasi-momentum K — ^|mnL| mo dL = and 
K = ~|(m + l)w£|modL = with m integer, respectively. Thus both of them 
belong to the K = block and are coupled by the interactions. On the other 
hand, in the incommensurate case, we have N — nL + AN, and hence the two 
states \N, 0) and |0, N) belong to different blocks and remain degenerate as K = 
\mnL + mAN\ modL ^ |(m + l)nL + (m + l)AN\ modL . 

In the commensurate case, we can construct an effective 2x2 Hamiltonian by 



projecting the many-body Hamiltonian (2.3) onto the subspace spanned by the 
states \N,0) and \0,N) 



H- 



E 



(i) 



2x2 



E 



(i) 

N 



where the coupling A is given in perturbation theory by 



A = 



E 



HoiHij 



pN 



p ( E Q 1) - £ i)( E Q 1 ' ~ £ j) ■ ■ ' ( E 0~' - e v) 



(1) 



u 

2L 



n(i-i) 



(2.5) 



(2.6) 



In this expression, n = N/L is the number density, the are transition matrix 
elements introd uced by the interaction Hamiltonian H{ nt , and e$ are either given 



by En^ in Eq. (2.4) or non-interacting many-body eigenenergies depending upon 



whether the intermediate states are or are not in the \n, N — n) manifold. The factor 
n(L — 1) corresponds to the minimum number of collision processes necessary to 
couple the states \N,0) and |0, JV) and the sum is taken over all possible coupling 
paths. In the case of the ring lattice the coupling is exclusively though states outside 
the degenerate manifold (see Fig. |3|, and the coupling A decreases exponentially 
with increasing number of particles N. 

At the critical phase twist A8 = it, w C nave — ±^ N 
zero value of the coupling A the symmetric and anti-symmetric superpositions |±) 
become the ground and first excited state separated by an energy gap 2A 



7r, we have E^ — EjP and due to the non- 



± 



| AT, 0) ± |0,iV) 

~7! ' 



(2.7) 



In Fig. [4] we plot the overlap with the states |iV, 0) and |0, A 7 ). Below and above 
the critical phase twist the bosons form the non-rotating condensate | N, 0) and the 
vortex state |0, N), respectively, while at resonance the macroscopic superposition 



state (2.7) occurs. These Schrodinger cat states are central to high-precision spec- 



troscopy, amplified quantum detection and measurement (Leibfried et al. 2004 



Cappellaro et al. 2005 Lee & Khitrin 2005 ) where they improve the resolution by 
a factor of VN with respect to the classical shot noise limit. 
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(c) Superposition states in rotating ring superlattices 



In Nunnenkamp et al. (2008) we further explored whether the situation can 
be improved by introducing a lattice modulation and considered instead a ring 
superlattices where Jj = J for j even and Jj — t for j odd. Due to the superlattice 
potential the quasi-momentum states \q) and \q + L/2) are coupled and the single- 
particle Hamiltonian is no longer diagonal in the quasi-momentum basis 

-(J + t) cos<j) —i(J — t)sm(j) \ f b q \ 



H*p- 2 ( % K+L2 J[ Jrl[J . /)sm ,, + ,./-:-/) cos,, )\ h t 

(2.8) 



with = 0- ^2. We can diagonalize the single-particle Hamiltonian rt2.8[) via a 
unitary basis transformation (c 9 ,c g+ ^/ 2 ) = Mu^,!),^) an d obtain 



where the single-particle energies are given by E^ = ± \J J 2 +t 2 + 2Jt cos (29 — ^jf) ■ 

In the uniform ring, t/J — 1, the eigenstates of the single-particle Hamilto- 
nian H sp are quasi-momentum states. At certain phase twists 9 they are doubly 
degenerate. For example, for L = 4 sites the quasi-momentum states \q — 1) and 
\q = —1) are degenerate at 6 = 0, whereas at 6> = ir/4 the states |g = 0) and 
\q = 1) as well as \q = 2) and \q = —1) are degenerate. Reducing the symmetry of 
the ring by choosing t ^ J the quasi-momentum states which differ by L/2 quasi- 



momentum units are coupled by the single-particle Hamiltonian (2.8), so that the 
quasi-momentum states \q — 1) and \q = — 1) hybridize and the degeneracy at 9 = 
is lifted, i.e. E+ - £f = 2( J - i). 

At 9 = 7r/4 however the degenerate quasi-momentum states are not coupled by 
the single-particle Hamiltonian ( |2.8| , so that the degeneracy is present also in the 
non-uniform case. This remains true for arbitrary L, i.e. Eq = E~^, 4 at 9 = 7r/4, 
but for L 7^ 4 these states are not the ground states of the system. In Fig. [T] we 
plot the single-particle spectrum for L = 4 sites as a function of the effective phase 
twist 9. It shows level crossings at 9 = ir/4 both for t = J as well ast^J. We will 
refer to 9 = 7r/4 as the critical phase twist, since - as we will demonstrate below - 
weak on-site interactions lift the degeneracy at 9 = ir/4 and lead to the formation 
of strongly-correlated states in the many-body system. 

The crossing of two single-particle levels implies a (iV + l)-fold degeneracy in the 
non-interacting many-body spectrum. In Fig. [2] we plot the many-body spectrum 
with L = N = 4 and t/J = 0.7 for U/J = and U/J = 0.5 as a function of the 
phase twist 9 and find that interactions turn level crossings into avoided crossings. 

In the weakly interacting regime this effect can be understood by deriving an 
effective Hamiltonian within the (iV+l)-dimensional degenerate subspace. A conve- 
nient basis for this subspace is spanned by the Fock states |n, N—n) with < n < N, 
where n particles are in the single-particle state of energy Eq and N — n particles in 
the one of energy res P ec tively. For weak interactions NU/L -C 2\J J 2 + t 2 this 

subspace is the low-energy sector of the many-body problem for all phase twists 9 



and tunneling strength ratios t/J. Starting from the interaction Hamiltonian (2.3) 
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we restrict the Hilbert space to the relevant modes and keep only terms within the 
low-energy subspace. In this way we obtain the effective Hamiltonian to first order 
in the on-site interaction strength U 



(e 



n 



U 



+ E L/4 n L/i )+— (2n n L/i + N 2 - N) + 



itjU 
~2L 



44 £ L/4CL/4 + H.C. 



where h 



c q c q are the number operators and the parameter rj evaluated at 9 



(2.10) 

tt/4 

is given by rj — jsxfi- The first bracket of Eq. (2.101 contains the contributions 
from the single-particle Hamiltonian (2.9), whereas the terms in the second and 
third brackets arise from the on-site interaction. At the critical phase twist 9 = ir/A 
the former are an unimportant zero-energy offset, whereas the terms in the second 
bracket shift the energies of the states in the subspace differently, e.g. they lead 
to an energy difference of U(N — 1)/L between the states \N,0) and \N — 1,1), 
while the states \n,N — n) and \N — n,n) remain pairwise degenerate. The terms 
in the third bracket are off-diagonal in the Fock basis of the subspace and describe 
two-particle scattering between the two single-particle modes. 

Let us now determine the ground and first excited state for slightly non-uniform 
rings t/Jfn 1, close to the critical phase twist 9 ps 7r/4. Since the terms in the sec- 
ond bracket of Eq. ( |2.10[ ) increase the energy for all states in the subspace apart 
from | TV, 0) and |0, N) and the coupling between the states is weak (as the coupling 



•q is small in this limit), we project the effective Hamiltonian (2.101 onto the sub 



space spanned by these two nearly-degenerate lowest-energy states. As there is no 
direct coupling between \N, 0) and |0, N) we calculate the total coupling through 
intermediate states using perturbation theory. After eliminating the intermediate 
states we obtain the following 2x2 Hamiltonian 



#2x2 — 



AE/2 
A* 



A 

-AE/2 



(2.11) 



where AE is the energy difference between the states \N,0) and |0, N) caused 
by the detuning of the phase twist from resonance A9 = 9 — i.e. AE = 
N(El /4 - Eq ) fa yffiffi , and A is the coupling between the states \N,0) and 
|0, N) due to the off-diagonal terms of the effective Hamiltonian (2.10). As the 
latter only directly couples the states \n, N — n) and \n ± 2, N — n =p 2), the first 
non-vanishing order is given by 



A 



(iV,0|i^ /2 |0,AO 



(1) 



E, 



u 

L 



N/2 



Nl 



(2.12) 



with the interaction energy shift E„ = |£ {2n(N — n) + N 2 — N) . Note that in 
contrast to the case of the ring lattice the perturbation couples states within the 
degenerate manifold (see Fig. [3]). This leads to a less severe but still exponential 



scaling of the gap A with the number of particles N (Nunnenkamp et at. 2008). 



The ground state of the two-by-two Hamiltonian (2.11) is similar to the one 



obtained above in Eq. (2.7), i.e 



a\N, 0) +i N / 2 (3\0,N) 
71 ' 



(2.13) 
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We see that to obtain a cat-like superposition, i.e. a//3 « 1, the energy difference 



must not dominate over the coupling \AE\ <c |2A| (see |Hallwood et al. (20071 and 



Nunnenkamp et al. (2008) for further details) 



3. Vortex nucleation in the lowest Landau level 



In this section, we discuss the effects of rotation on ultracold bosons confined to a 



two-dimensional harmonic potential. We hrst review the results of Dagnino et al. 
(20096 ah in order to then compare and contrast them with our findings on rotating 



ring (super)lattices which we presented in the previous section. 



(a) Hamiltonian 



Following Dagnino "eFaT1p0096T 



we consider a system of N ultracold bosons 
with mass M confined to a two-dimensional symmetric harmonic potential Vq and 
rotating in the x — y plane about the z-axis with angular velocity 51, then as before 
the many-body Hamiltonian in the rotating frame is 



H = / dx& 



2M 



v 2 + y (x) + 



Airh 2 a s 
2M 



(3.1) 



If we assume that the potential Vq confines the motion along the z-axis so 
strongly that only the transversal motion of the atoms in the x — y plane needs to 
be considered, and that along the transverse direction the atoms feel a harmonic 
confinement Vo(x,y) = Muj 2 (x 2 + y 2 )/2 and the interaction energy is small com- 
pared to the Landau level splitting h(uj + f2), the bosonic field operator <& can be 
expanded in the lowest Landau level basis $ = ^ m a m (p m (x,y). Here, ip m (x,y) 
are the eigenfunctions of the single-particle angular momentum operator L z with 
non-negative integer eigenvalue m, i.e. <p m (x,y) cx (x + iy) m e~^ x +v ^ 2X with the 
magnetic length A = \JK[Moj. In the course of the following discussion an asym- 
metry of the trapping potential in the x — y plane will be included by adding a 
single-particle potential V(x,y) = 2AMuj 2 (x 2 — y 2 ), where A is a measure of the 
asymmetry. For A ^ 1 the lowest Landau level remains a good basis set and the 
asymmetry can treated as a perturbation. Within this approximation the many- 

2009fe|c 



body Hamiltonian is (Dagnino et al. 



H = H + U + V 



h(u - Q)L + 



A 



E 



(k + l)\S i+ j k+i ,+-t- - 
a- a^a k ai 



4ttA 2 ^ 2 k + l Ji\j\k\l\ 

ijkl 



+ — A 2 ^ \/ m ( m - l)aj n a m _ 2 + \]{m + l)(m + 2)a) m a m+2 . 



(3.2) 



Here, Hq is the unperturbed single-particle Hamiltonian proportional to the total 
angular momentum operator L = J2m ma m a m, U is the two-body interaction and 
V is the perturbation due to the asymmetry in the single-particle potential. In the 
following we use A, Huj and oj as units of length, energy and frequency, respectively. 
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(b) Effective many-body Hamiltonian 

Analogous to the previous section, we start our discussion with the single- 
particle spectrum E m = h(uj — fl)m. In Fig. [I] we plot it as a function of rotation 
frequency f2. We see that all states in the lowest Landau level become degenerate 
at Q — u>, but in contrast to the rotating ring lattice there is no single-particle level 
crossing for < £1 < w. 

Nonetheless, the spectrum of N non-interacting bosons is highly degenerate for 
any rotation frequency il. This is a consequence of the fact that the single-particle 
Hamiltonian is not only diagonal in the angular momentum basis but proportional 
to the total angular momentum L, so that all many-body states with the same 
total angular momentum L have the same energy - a degeneracy independent of the 
rotation rate Q and only present at the many-body level. Interactions will break this 
huge degeneracy, but since the interaction Hamiltonian U commutes with the total 
angular momentum operator L the interacting many-body Hamiltonian remains 
block diagonal, i.e. it only mixes states with the same total angular momentum L. 



Smith & Wilkin (20001 have shown that the states with the wave functions 



®L(xi,yi,...,x N ,y N ) oc (z i± - z ) . . . (z tL - z ) $ (x 1 ,y 1 , . . . ,x N ,y N ) 

l<i ± <...i N <N 

(3.3) 

where Zj — Xj + iyj are complex coordinates in the x — y plane, zq — y\ Zj/N is 

the center of mass, and $ oc e i' , i.e. the non-rotating ground-state wave 
function, are exact eigenstates of the interaction Hamiltonian U with eigenvalue 

E L = y — {2N - L - 2) . (3.4) 

They also present overwhelming numerical evidence that these are the ground states 
for < L < N with L ^ 1. When expressed in the Fock basis \no,ni, . . .) with 
occupation numbers n m for the angular-momentum basis, these are seen to be 
complicated superposition states. For example, the many-body state with one unit 
of angular momentum per particle |$iv) is not, in fact, the vortex state |0, N, 0, . . .) 
with all particles occupying the single-particle wave function with m = 1, but rather 
the so-called yrast state whose occupation number of the first angular-momentum 
state is large (ni) ~ N but whose overlap with the vortex state is only about one 



half fBertsch & Papenbrock| 


1999). 


Following 


Dagnino et al. 


20096 



a function of rotation frequency Q. Since the single-particle and interaction energy 
of the states \3>l) depend linearly on L, there is a rotation frequency fl c /u = 1 — 
at which all many-body states are degenerate. Around this rotation frequency 
the ground state rapidly changes: for fl < £l c the ground state is |$o)j while above 
SI > S7 C it is the yrast state |$iy). 

Turning to the effect of the perturbation V on the states involved, we plot in 
Fig. [T] the single-particle spectrum as a function of rotation frequency f2. We see 
that apart from a single-particle crossing close to Q w the perturbation has little 
effect on the single-particle spectrum, and it turns out that this level crossing is 
unimportant for the discussion that follows. 

The many-body spectrum, on the other hand, is degenerate at Q c , so that any 
additional perturbation can effectively mix the degenerate many-body states 
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Like Dagnino et al. (2009&), we show in Fig. [2] the many-body spectrum for A =/= 0, 



where we find the degeneracy is lifted into a set of anti-crossings 

At first glance the situation seems to be similar to the case of the superlattice 
discussed in the previous section. In both cases the N + 1-fold (iV-fold) degeneracy 
in a set of many-body states is lifted by an external perturbation (interactions 
in the superlattice case and trap asymmetry in the LLL case) which couples the 
various many-body states among them. This is not the case in the uniform ring 
lattice where coupling takes place via non-resonant states (see Fig. [3]). 

However, there are two crucial differences between the two systems. First of all, 
the interaction Hamiltonian U does not simply couple the degenerate states but also 
shifts their energies differentially (see Fig. [3 1. Th is is why we were able to construct 
an effective two-by-two Hamiltonian i?2x2 (2.11 1 that only couples the two lowest- 
lying states, and found the ground-state to be a macroscopic superposition state 
exhibiting an energy gap which exponentially decreases with increasing number 
of particles (due to virtual couplings to the adiabatically eliminated states within 
the resonant manifold). In contrast, the perturbation V contains only off-diagonal 
matrix elements and all states within the degenerate manifold have to be treated on 
an equal footing. This view is corroborated in Fig. [4] where we, like [Da gnino et al. 



(2009a), plot the overlap of the ground state with the states \<&l) as a function of 
rotation frequency fi. We see that with increasing Q the weight of the ground state 
in the various states shifts towards higher L in steps of two units of angular 



momentum. At the nucleation point Dagnino et al. (2009&) find that the state of 



the system is well described by a complicated superposition state 

|*o> oc \N, 0) + \N - 2, 2) + • • • + \2, N - 2) + |0, N) (3.5) 

where \n, m) is the state with n and m atoms in the two eigenfunctions of the single- 
particle density matrix wy>(r,r') = {^q^S!^ (r)^{r' with the largest eigenval- 
ues, respectively. The nature of the intermediate states also leads to a different 
scaling of the energy gap with the number of particles. As is shown by [Dagnino 



et al. (20096 1, the energy gap within the subspace at the critical rotation frequency 
remains finite as the number of particles increases, as opposed to the exponential 
scaling of the energy gap for the macroscopic superposition states in the ring lat- 
tice. Note, however, that the energy difference between the ground state within the 
subspaces of even and odd L decreases exponentially with the number of particles 



(Parke et al. 2008) 



The second difference between vortex nucleation in the two systems is the ground 
state above the critical rotation frequency. In the LLL system it has a large overlap 
with the yrast state | jv) - T° compare with the superlattice system we plot in 
Fig. [4] also the overlap between the ground state and the macroscopic occupied 
modes on both sides of the resonance: i.e. the condensate at rest \N, 0, ...) for 

< ri c and the vortex state |0, N, ...) for £1 > fl c . We see that whereas the 
overlap with the corresponding states is always close to one for the superlattice, in 
the LLL the overlap is only close to one for f2 < f2 c . For fi > fl c is only about one 
half. This is a consequence of the fact that in the ring lattice away from the critical 
rotation frequency the state |0,./V) is the only low-lying state and therefore finite 
interactions lead to a small depletion of the superfluid ground state. In contrast, 
the non-interacting LLL system with one unit per particle is highly degenerate and 
interactions produce strongly-correlated many-body ground-states like | jv) • 
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4. Noise Correlations 



To further quantify the entanglement of the ground state at the critical frequency 
and to compare the different nature of the ground state in the ring lattice and LLL 
system we calculate the noise correlation pattern, i.e., 



(d\,6 q ,)(6\6 q ). 



(4.1) 



We will use quasi-momentum creation operators for the lattice case O q = b q and 
angular- momentum operators for the harmonically trapped system, O q = a q . 

(2009) have argued that noise correlations are experimentally 



Mintcrt et al. 



accessible quantities encoding information on entanglement. An alternative char- 
acterization of the entanglement properties of the ground state in the LLL system 



has been given by Liu et al. (2009), and Read & Cooper (2003) have proposed 



time-of-flight expansion to probe the vortex lattice in LLL systems. 

For the lattice and superlattice cases the noise interferogram shows the de- 
velopment of three sharp fringes, with positive and negative signs at the critical 
frequency. Their amplitudes are given by 



A(0, 0) = A(l, 1) = -A(0, 1) = N 2 /2 



(4.2) 



heralding the development of a Schrodinger cat state at Q, c . Away from f2 c the noise 
interference pattern disappears signaling the unentangled nature of the ground state 



and the macroscopic occupation of a single mode ( |Rey et al. 2007 ) 



In striking contrast is the noise interferogram for the LLL system which exhibits 
various sharp fringes which develop as the system is driven through the critical 
frequency. They signal the strong correlations and multi-mode nature of the ground 
state at fi c . Furthermore, the interferogram does not disappear for f2 > fl c , behavior 
which highlights the correlated nature of the yrast state |<J>at). All these features 
are illustrated in Fig. [5] where we plot the non-zero A(q, q') as a function of phase 
twist 9 and rotation frequency f2, respectively. 



5. Conclusion 

We have shown how a new generation of experiments can be used to examine the 
quantum nucleation process and its important link to entangled states of atoms. 
We have compared quantum vortex nucleation in rotating ring lattices and two- 
dimenisonal harmonic potentials and have shown how effective Hamiltonians can be 
used to describe the nearly-degenerate ground-state manifold close to the rotation 
frequency at which the first vortex is nucleated. The degeneracy in the many-body 
spectrum is lifted by interactions in the case of the ring superlattice and by the 
asymmetry of the single-particle potential in the LLL system. In the first case the 
interactions not only couple the states but also shift their energies differentially. This 
leads to macroscopic superpositions as low-lying states whose energy gap decreases 
exponentially with increasing number of particles. In contrast, in latter case all 
states contribute to the ground-state wave function and the energy gap within the 
subspace remains finite with increasing particle number. Finally, we showed that 
the two scenarios can be distinguished in noise correlation interferograms: while 
the ground state in a rotating ring lattice above the critical phase twist is the 
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Figure 1. Single-particle spectrum versus phase twist (or rotation frequency). 

(Upper) Rotating ring lattice with L — 4 sites: there is a single-particle level crossing at 
6 = tt/4 both for t = J (left) and t = 0.7 J (right) (Nunnenkamp et al.| [2008[ ). (Lower) 
LLL system: both for A — (left) and A = 0.03 (right) the single-particle spectrum is 
degenerate only at fi = U) for A = (left) or at SI fa uj for A — 0.03 (right). This crossing 
at single-particle level is not important for the discussion that follows. The insets show 
the details around fi « w. 



vortex state with vanishing noise correlations, interactions in the lowest Landau 
level produce the yrast state which has a non-trivial noise pattern. 
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Figure 3. Coupling path. The degeneracies of the many-particle spectrum are lifted by 
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tively) or by the asymmetric potential V in case of the LLL system (lower). While in the 
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Figure 4. Properties of the ground-state wave function. (Left) Overlap of the ring 
lattice ground state with the wave function of a condensate in the two lowest-lying quasi- 
momentum states: away from the critical phase twist the atoms condense (Rey et al 



20071. (Right) Overlap of the LLL ground state for N = 6 particles in the asymmetric 
case (g — 1 and A — 0.03) with the basis set of the effective Hamiltonian: it changes 
from |$o) via the superposition (3.5 I to the yrast state |$jv) (Dagnino et al. 200961. We 



use exact diagonalization of the many-body Hamiltonian (3.2 1 including all states up to 
a total angular momentum of L max — N + 2. For comparison the inset shows the overlap 
of the LLL ground state with non-rotating condensate (solid) and vortex (dashed) states. 
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form a strongly-correlated many-body state. 
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